Partition Function Zeros of a Restricted Potts Model on Self-Dual Strips of the 

Square Lattice 
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We calculate the partition function Z(G,Q,v) of the Q-state Potts model exactly for self-dual 
cyclic square-lattice strips of various widths L y and arbitrarily great lengths L x , with Q and v 
restricted to satisfy the relation Q — v 2 . From these calculations, in the limit L x — *• oo, we determine 
the continuous accumulation locus B of the partition function zeros in the v and Q planes. A number 
of interesting features of this locus are discussed and a conjecture is given for properties applicable 
for arbitrarily great width. Relations with the loci B for general Q and v are analyzed. 
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I. INTRODUCTION 

The Q-state Potts model has served as a valuable 
model for the study of phase transitions and critical phe- 
nomena On a lattice or, more generally, on a 
graph G with n sites (vertices), at temperature T, this 
model is defined by the partition function 



Z(G,Q,v) = ^exp(if^<5 aiaj ) 



(1.1) 



M 



where cr, = 1, q are the classical spin variables on each 
vertex i S G, (ij) denotes pairs of adjacent vertices, K = 
PJ where j3 = (/c^T) -1 , and J is the spin-spin coupling. 
We define 
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(1.2) 



so that v has the physical range of values < v < oo and 
— 1 < v < for the ferromagnetic and antiferromagnetic 
cases J > and J < 0, respectively. The (reduced) free 
energy per site is / = — (3F = linin-^oo n" 1 In Z(G, Q, v). 
Aside from the Q — 2 Ising special case on two- 
dimensional lattices, the free energy of the model has 
never been computed exactly for general Q and v on lat- 
tices of dimensionality d > 2, although the critical behav- 
ior is understood for d — 2. In view of this lack of exact 
solutions for d > 2, there is continuing value in studying 
the case of lattice strips of finite width L y and arbitrarily 
great length L Xl with various boundary conditions. Al- 
though these are quasi-one-dimensional systems, and the 
resultant free energy does not exhibit a phase transition 
at any finite temperature, they are nevertheless useful for 
the exact results that can be obtained and the insights 
that one can gain concerning the properties of the model 
as a function of Q and v. Among these exact results are 
the partition function zeros and the loci that they form 
in the limit n — > oo. These loci are determined by the 
equality in magnitude of the eigenvalues of the transfer 



matrix with maximal modulus and hence are also called 
the set of equimodular curves. (Here "curves" is used is a 
generic sense; these loci can also include line segments.) 
In previous works we have presented exact calculations 
of the Potts model partition function for arbitrary Q and 
v on such strips and have also investigated various spe- 
cial cases 0. One type of special case is provided by 
considering requiring Q and v to satisfy some functional 
relation. 

In this paper we consider the Potts model on self-dual 
cyclic square-lattice strips of various widths L y and ar- 
bitrarily great lengths L x , with the variables Q and v re- 
stricted to satisfy, in terms of the variable Q q = v/\/Q, 
the relation 
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sq 
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(1.3) 



It will actually be convenient to consider the square of 
this condition, namely, 



Q = v 2 



(1.4) 



The condition l|1.3fl is the condition for the phase tran- 
sition temperature of the Potts ferromagnet on a two- 
dimensional square lattice, i.e., kgT c = J/ln(l + ^/Q ). 
(The transition is second-order for < Q < 4 and first- 
order for Q > 4.) For the L x — > oo limits of these strips 
we exactly determine the continuous accumulation loci B 
of the partition function zeros in the v and Q planes. 

Because of the quasi-one-dimensional nature of the 
infinite-length strips that we consider, the significance 
of the condition (|1.3(l or 1)1. 4JI is not that the Potts ferro- 
magnet has a phase transition at this temperature, but 
rather that it is a self-dual point; the partition function 
satisfies the duality relation 



Z(G,Q,v) oc Z(G*,Q,v d ) , 
where G* denotes the planar dual to G and 
Q . 1 



i.e. 



C 



(1.5) 



(1.6) 



If G = G* , as is the case with our self-dual square-lattice 
strips, and if v = v^, then, up to a multiplicative factor 
which will not be important here, the partition function 
is transformed to itself under duality. The properties of 
the Potts model for the submanifold of values of Q and v 
that satisfy the condition (|1 .4|> are thus of general interest 
even if the model does not exhibit a phase transition at 
the corresponding temperature. 

II. PRELIMINARIES 

In this section we review some background which will 
be used in the subsequent portion of the text. A graph 
G = G(V,E) is defined by its vertex set V and its 
edge set E (where V should not be confused with the 
temperature-like Boltzmann variable v). The number n 
of vertices of G is denoted as n = n(G) = \V\ and the 
number of bonds (edges) of G as e(G) = \E\. The de- 
gree of a vertex is defined as the number of edges that 
connect to it. The Potts model can be generalized from 
non-negative integer Q to non-negative real Q via the 
relation Q 

Z(G,Q,v) = J] Qfc(G') t; e(G') (2.1) 
G'CG 

where G' is a spanning subgraph of G, i.e., G' = (V,E') 
with E' C E, and k(G') denotes the number of con- 
nected components of G' . This relation also allows one 
to consider v and Q to be complex, as will be necessary 
in analyzing the zeros of Z(G,Q,v). A transfer matrix 
formalism for the Potts model that can be used for real 
Q was presented in Ref. @|. 

Given that one restricts Q and v to satisfy eq. 
one can express Z{G, Q, v) as a function of one variable. 
Expressing Q as a function of v via eq. eq. I|1.4(l . one 
obtains Z(G, v 2 ,v). This has the general form, from eq. 

en, 

Z(G,v 2 ,v) = J2 « 2fe(G ' )+e(G,) (2.2) 

G'CG 

The alternate procedure, expressing v as a function of 
Q, is not as convenient, since the formal solution for 
v involves a sign ambiguity, v = ±y/Q~. Clearly, since 
Z(G, v 2 , v) is a polynomial in v with real coefficients [jj, 
the set of its zeros and the resultant locus B are invariant 
under complex conjugation. 

The Potts model partition function is equivalent to an 
important function in mathematical graph theory, the 
Tutte polynomial T(G,x,y) HQ: 

Z(G, Q,v) = (x- l) k ^ (y - 1)"(°) T(G, x, y) (2.3) 
where 

x = 1 + — , y = 1 + v . (2.4) 



(We follow standard mathematical notation here; the 
reader should not confuse these Tutte variables x and 
y with the longitudinal and transverse strip graph coor- 
dinates x and y, since the meaning will be clear from 
context.) For a planar graph G, the Tutte polynomial 
satisfies T{G, x, y) — T(G* , y, x). In terms of Tutte vari- 
able, eq. (|1.4|) reads (x — l)(y — 1) = (y — l) 2 . Since the 
case y = 1, i.e., K — 0, is trivial, one may concentrate 
on the case y ^ 1 and divide by (y — 1), obtaining the 
corresponding condition x — y. For a self-dual planar 
graph G such as the lattice graphs that we consider here, 
T(G,x,y)=T(G,y,x). 

The Potts model on a L x x L y section of a lattice is 
related to quantum algebras, in particular, to the algebra 
defined by a set of operators Ui satisfying the relations 

m 

U 2 = Q l/2 U l (2.5) 
[U i ,U j ] = Q if (2.6) 

UiU i±1 Ui = Ui (2.7) 

where the range of the index i on Ui can be taken to be 
1 < i < 2L y — 1. One commonly writes 

VQ = q + q^ (2.8) 

where 

q = e ie < 2 (2.9) 

so that 

Q = 4cos 2 (6»/2) . (2.10) 

The case where q is a root of unity, i.e., 9 /it is rational, 
will be relevant to our discussion below. We define the 
root of unity q r = e m l r and correspondingly 

Q r = 4cos 2 (7r/r) . (2.11) 

The values Q r corresponding to the roots of unity q r 
are often called Tutte-Beraha numbers; early interest 
in these arose in studies of the chromatic polynomial 
P(G, Q) = Z(G, Q, —1), equal to the Potts antiferromag- 
net at zero temperature [I HE1 (see also P2). They 
also play a special role for the finite-temperature Potts 
model (with either sign of J), as is clear from the fact 
that the representations of the algebra of eqs. I|2.5|l - I|2.7|l 
are reducible at Q = Q r 00. We also define 

v r = -2cos(7r/r) (2.12) 

so that Q r = v 2 . For Q = Q r and v = v r , the associ- 
ated values of the Tutte variables are given by eq. Ij2.4|l 

as Xr = yr 

= 1 — 2cos(7r/r). Calculations of zeros of 
Z(G,v 2 ,v) on finite sections of the square lattice found 
that a subset of these zeros occurred at or near certain 
Q r values an d, as will be discussed below, we find 
that the accumulation loci Bq cross the real Q axis at 
certain Q r values. 
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III. GENERAL FORM OF Z(G,Q,v) FOR 
SELF-DUAL CYCLIC SQUARE-LATTICE STRIPS 

We focus on strip graphs of the square lattice of length 
L x vertices and width L y vertices with periodic longitudi- 
nal boundary conditions and with one additional external 
vertex connected to each of the vertices on one side of the 
strip, say the top side. For these graphs the number of 
vertices is n — L x L y + 1, which is equal to the number 
of faces, and the number of edges is e = 2n 15] . The 
average vertex degree, defined as 2e(G)/n(G), is equal to 
4, the same value for the square lattice. As regards indi- 
vidual vertices, the external vertex has degree L x , the L x 
vertices on the bottom side of the strip have degree 3, and 
the L x (L y — l) other vertices have degree 4. These graphs 
are self-dual, and accordingly we label the boundary con- 
ditions as self-dual cyclic boundary conditions (CDBC) 
and the graphs as Gcdbc[L v x L x \. The fact that these 
square-lattice strip graphs maintain the property of self- 
duality that holds for the (infinite) square lattice makes 
them particularly valuable to use in this type of study. 
Furthermore, the periodic longitudinal boundary condi- 
tions minimize boundary effects in this direction. To ex- 
plain one of the most important consequences of these 
self-dual boundary conditions, we recall that the locus 
Bq for the Potts ferromagnet on infinite-length, finite- 
width strips usually extends to £ = oo, i.e., K = oo, 
reflecting the fact that the model has a critical point at 
T = (and zero magnetization for T > but nonzero 
magnetization at T = 0). This was, indeed, the case for 
strips of many different lattices for which we obtained ex- 
act determinations of the locus Bq or the equivalent loci 
B v or B a (e.g., 24]- 29]). However, as our exact determi- 
nations of Bq for infinite-length, finite- width strips of the 
square lattice with cyclic or free self-dual boundary con- 
ditions in Ref. 21] showed, the locus B^ does not extend 
to £ = oo and is compact in the £ plane. This difference 
arises because, owing to the self-duality, it follows that 



form \M\1 



Bq is invariant under £ — > — 



Consequently, for the self-dual strips, B^ passes through 
the zero-temperature point £ = oo if and only if it passes 
through the infinite-temperature point, £ = 0. But it 
cannot pass through the latter point because the free en- 
ergy is analytic in the neighborhood of K = 0. Thus, for 
general Q and v, the infinite- length, finite- width self-dual 
strips share with the (infinite) square lattice both the 
property of self-duality and the property that B^ is com- 
pact in the £ plane. As we shall show, however, when we 
restrict Q and v to obey the relation l|1.4|) . the resultant 
singular locus B v , and hence also Bq, are noncompact 
in their respective v and Q planes. We note that finite 
sections of the square lattice with these CDBC bound- 
ary conditions have previously been used for studies of 
complex-temperature zeros of the Potts model partition 
function in Refs. [Il[Il[l9|. 

The Potts model partition function has the general 



Z{Gcdbc[L v x L x ],Q,v) 
Qv 



Ly + l n z {L y ,d) 

L:Lv E E S^fo,,^)*" (3.2) 

d=l j=l 



where 



K w^Q- K w = Y j{ -iy( 2d - 1 - i \ Q d- 

j=0 v J 



and 



n z {L yi d) 



2d ( 2L y + 1 



L y + d+l \L y - d + 1 



(3.3) 



(3.4) 



for d = 1, L y + 1 (and zero for other values of d). The 
^L y ,d,j are eigenvalues of a transfer matrix of the type 
discussed in Ref. \&. In terms of the notation of Ref. 

V = '• ' A,. • ( 3 - 5 ) 

Since nz(L y ,L y + 1) = 1, i.e., there is only one eigen- 
value for d = L y + 1, we label it simply as Xl l v +i,i = 
Al„,l b +i. This eigenvalue is [MEH 



(3.6) 



From the general form l|3.2|l it follows that 
Z(GcDBc[L y x L x ],v 2 ,v) has a zero of multiplic- 
ity at least L x L y + 2 = n+ latw = 0. We find that this 
is an isolated zero for these types of graphs. 
The coefficients can be written as 120, 23 1 



= l[(Q-s d , k ) 



(3.1) where 



fc=i 



s d ,k = 4 cos 2 ( ^ ) 



2d J 



(3.7) 



(3.8) 



The first few coefficients «w are = 1, k' 2 ' = Q — 2, 
and = (Q - 1){Q - 3). From eq. ESJ or ^ it 
follows that if d > 3 and d = mod 3, contains the 
factor (Q — 1)(Q — 3). i.e., 

R {d) = {Q - 1){Q - 3) fi(Q) if d>3, d = mod 3 

(3.9) 

where [i(Q) is a polynomial of degree d—3. This property 
will be used later. We also note that if d is even, then 
/«w contains the factor [Q — 2). The total number of 
eigenvalues is 



J2n z (L y ,d)-(" 2L -' ' ' 
d=i 



Ly + l 



(3.10) 
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For given values of Q, v and d, the eigenvalue with maxi- 
mum modulus, i.e., the maximal or dominant eigenvalue, 
is labelled \ Ly ,d, ma x- 

As discussed in (section 2.2 of) Ref. [24|, the formal 
definition for the reduced free energy, / is not adequate 
for (real, non-negative) non-integer Q because at certain 
special points Q sp one has the noncommutativity of limits 

lim lim Z{G,Q,v) 1/n ^ lim lim Z{G,Q,v) 1/n . 

n— >oo Q-*Q 3p Q^Qsp n—>ao 

(3.11) 

The definitions of / with these two different order of lim- 
its are denoted, as in Ref. |24j|. f n g and /q„. Clearly, 
no such issue of noncommutativity arises if one restricts 
to positive integer Q values and uses the original Potts 
model definition l|l.l|l . However, since we will consider 
the more general case of positive real Q as encompassed 
in the formula 1)2.1(1 . we will have to take account of this 
noncommutativity. 



IV. PARTITION FUNCTION ZEROS AND LOCI 

B 

A. General 

We now determine the loci B in the v and Q planes, 
denoted B v and Bq, for the infinite- length limits of these 
strip graphs for several widths L y . We also compare these 
to zeros calculated for long finite strips. Our procedure 
is to calculate Z(Gcdbc[L v x L x ],v 2 ,v), solve for the 
locus B v as the locus of solutions of degeneracy in mag- 
nitude of dominant Ai dj's, and then obtain Bq via the 
mapping ((1.4J1 . We find that these loci exhibit a number 
of interesting systematic properties. To explain these, 
we first discuss our exact solutions for strips of specific 
widths. 



In this case, the graph is equivalent to the wheel graph 
in mathematical graph theory. Using our exact solution 
for the Potts model partition function on this graph [2^, 
El El, we find that, with L x = m, 



Z(G C dbc[1 x m],v 2 ,v) 



„m+2 



{(A 1 , 1 , 1 r + (A 1 , 1 , 2 r} + s( 2 ) 



where 



Ai 



3 + 2v ± V5 + 4v 



(4.1) 



(4.2) 



with j — 1,2 corresponding to ±, and we have used 
A12 = 1 from cq. (|3.6[) . Here and below the coefficient 
functions = K^ d '(Q) are understood to be evaluated 
on the manifold (|1.4|) . 



In the limit L x — > oo, we have calculated the exact 
locus B v in the v plane and its image under the mapping 
(|1.4J) . Bq in the Q plane. We find that B v consists of the 
union of (i) a pinched oval that crosses the real v axis at 
v = U3 = — 1 and v = — 2 with (ii) a semi-infinite line 
segment on the negative real v axis extending leftward 
from v = — 2 to ?) = - oo. In the region outside of the 
pinched oval, is the dominant eigenvalue. Along 

the semi-infinite line segment, = XI 1 2 , so that 

|Ai,i.i| = | Ax . i ,2 1 - In the region within the oval, A1.2 = 1 
is the dominant eigenvalue. The fact that B v does not 
intersect the positive real axis at any finite point is a 
general result reflecting the quasi-one-dimensional nature 
of these infinite-length strip graphs and the consequent 
property that the free energy is analytic for all finite tem- 
peratures. As one moves leftward along the real v axis, 
the first nonanalyticity that one encounters, at v = — 1 
is precisely the zero-temperature limit of the Potts anti- 
ferromagnet (with Q = v 2 = 1). 

The locus Bq, shown in Fig. ^ is the union of (i) a 
pinched oval that crosses the real Q axis at Q = Q3 = 1 
and Q = Q\ =4 (where Q r was defined in cq. (12.11(1 ) 
with (ii) a semi-infinite line segment on the real Q axis 
extending from Q = 4 to Q = 00. Note that in terms of 
the Tutte variables, 



-2, Q 



y 



1 



(4.3) 



In Fig. Q]we also show partition function zeros calculated 
for a long finite strip, with L x — 40. (For this length some 
of the zeros on the positive real have magnitudes too large 
to be included in the figure.) At v — —2 and its image 
point Q = 4, the magnitudes of all three eigenvalues are 
degenerate and equal to 1. In addition to the isolated 
zero at Q = Q2 — as discussed above, the partition 
function has a zero very close to Q = Q4 = 2. This 
can be understood as a consequence of the fact that in 
this region X12 max is dominant, but its coefficient — 
Q — 2 vanishes at Q = 2. In Fig. n° ne can a l so see nas 
several real zeros in the interval Q > 4; as L x — > 00, 
these and other real zeros in this interval merge to form 
the semi-infinite line segment on Bq. 

This calculation provides a simple example of the non- 
commutativity (|3.11|) . Consider the point v sp = — v2, 
whence Q sp = v 2 p = 2. Now if we first set Q — 2 and 
then take n — * 00 to calculate f n Q, then, since the co- 
efficient vanishes at this point, the dominant term 
A" l 2 = 1 in eq. 14. 1|) is annihilated. Furthermore, at this 
point 



Ai. 



1,3 



l)e 



± i0 



where 



so that 



arctan =- 

V 3-2V^ 



(4.4) 



(4.5) 



Z(G C DBc[lxm},2, -V2) - 2(-2)^ (V2-l) m cos(mfl) 

(4.6) 
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FIG. 1: Locus Bq for the Potts model on a 1 X oo graph with CDBC 
and Q = v 2 . This locus includes a semi-infinite line segment on 
the interval 4 < Q < oo. Partition function zeros are shown for a 
graph of length L x = 40. 



and hence 

lim lim _ \Z 1 ' n \ = 2 - V2 (4.7) 

n^oo (Q jt) )_,(2,-V2) 

With the opposite order of operations, one obtains the 
different result 

lim _ lim \Z 1/n \ = V2 (4.8) 

(<3,t0->(2,— v/2) n ^°° 

where, since there is no canonical phase to pick for the 
1/n'th root (24|, we just display the modulus. 



For the width L y — 2 strips with CDBC we use our 
previous exact calculation of the partition function in 
Refs. [53, and specialize to the submanifold of vari- 
ables given by (|1.4|l . From eq. 1)3. 4f> there are five A2,i,j's, 
and four A2,2,/s, so that, together with the single A2,3, 
the total number of \2,d,j' s lS ten. We have 



r 5 



Z(G C dbc[2 xm],v 2 ,v)=v 



2m+2 



3=1 



where 



(4.9) 



1/2 



7+r ? \/5+4u± 2(19+7?yV5)+8(3+r;V / 5)t; 

(4.10) 



where r] — ±1, yielding four possible terms. The A2,i,j 
are roots of a quintic equation which can be derived from 
the general quintic given in Ref. [2l|. For v = —2, all 
of the roots (|4.1U|I have unit modulus and this is also 
true of the roots of the above-mentioned quintic equa- 
tion, so that all ten \2,dj have |A2,dj| = 1- We show 
our results for B v and Bq in Figs. El and together 
with partition function zeros calculated for a strip with 
length L x — 30. (For this and other figures, some zeros 
have sufficiently large magnitudes so that they are not 
shown.) We find that the locus B v crosses the negative v 
plane at v = -2 and v = -2 cos(-7r/(2£ + 1)) for £=1,2, 
i.e., -1, and -2cos(vr/5) = -(1 + Vb)/2 ~ -1.618. 
Correspondingly, Bq crosses the positive Q axis at the 
squares of these values, 4cos 2 (tt/(2£ + 1)) for I = 0, 1, 2, 
i.e., 4, 1, and 4cos 2 (vr/5) = (3 + v / 5)/2 ~ 2.618. Two 
complex-conjugate curves on B extend to complex infin- 
ity in the v and Q planes, rendering the loci noncompact 
in these variables. In contrast, B is compact in the 1/v 
or 1/Q planes, and we show B in the 1/Q plane in Fig. 0] 
The fact that B extends to 1/Q = means that the free 
energy does not have a Taylor series expansion in 1/Q. 
This feature is similar to the situation for the quantity 
W = linin^oo P(G, Q) x l n for certain families of graphs 
[3fj|. The confluence of curves at v — — 2 and thus Q = 4 
reflects the degeneracy in magnitude of the ten A2,d,/s. 

We comment on some additional features of Bq in 
Fig. 01 This locus contains a set of complex-conjugate 
arcs that extend from Q = 4 outward from the real 
axis, crossing the imaginary axis at Q ~ ±6.52i and 
ending at Q ~ —1.96 ± 6.27i. Along the real Q axis 
for Q > 4, the dominant term is A2,i jm aa;- A different 
A2,i, maa; is dominant in the region extending to complex 
infinity in the Q plane separated by curves on Bq from 
this wedge containing the interval Q > 4. In the re- 
gion containing the real interval Q3 < Q < Q5, i.e., 
1 < Q < (l/2)(3 + v/5) ~ 2.618, X 2 ,2,ma X is dom- 
inant, while in the region containing the real interval 
Q5 < Q < Qi, i-e., (l/2)(3 + \/5) < Q < 4, A 2 , 3 = 1 
is dominant. The partition function also has several iso- 
lated zeros that are not on Bq. In addition to the zero 
at Q = 0, it has zeros in the interval [0, 4] very near to 
Q = Qi = 2 and Q = Qq = 3. The zeros very near to Q4 
and Q e result from the fact that the coefficient that 
multiplies the dominant term (A2, 2. max)" 1 in the region 
around Q = 2 vanishes at Q = 2, and the coefficient 
that multiplies the dominant term ^2,3)™ = 1 in the 
region around Q — 3 vanishes at Q — 3. The partition 
function also has real zeros in the region Q > 4. 



D. 



For L y = 3, from the exact calculation of 
Z(Gcdbc[3 x L x ], Q,v) performed for Ref. [2l| we spe- 
cialize to the manifold (|1.4|l and calculate the partition 
function zeros. Because of the large number of A3 d 3 -'s (a 
total of 35) that enter, we do not list these explicitly here. 
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FIG. 2: Locus B v for the Potts model on a 2 X oo graph with 
CDBC and Q = v 2 . Partition function zeros are shown for a graph 
of length L x = 30. 



FIG. 4: Locus B plotted in the 1/Q plane, for the Potts model on 
a 2 X oo graph with CDBC and Q = v 2 . Partition function zeros 
are shown for a graph of length L x = 30. 
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FIG. 3: Locus Hq for the Potts model on a 2 X oo graph with 
CDBC and Q = v 2 . Partition function zeros are shown for a graph 
of length L x = 30. 



Having discussed above how the zeros in the v plane and 
the locus B v map to zeros and the locus Bq in the Q 
plane, we display our results only in the latter plane. In 
Fig. [SI we show a plot of zeros and Bq. To depict the por- 
tions of Bq that extend infinitely far from the origin, it is 
again useful to include, as Fig. a plot of B in the 1/Q 
plane, where it is compact. We find that in the interval 
— 2 < v < — 1 the locus B v crosses the negative real v 



axis at v = -2 and v = -2 cos(tt/(2£ + 1)) for £ = 1,2,3, 
i.e., -1, -2cos(tt/5) = — (1/2)(1 + y/E) ~ -1.618, and 
— 2cos(7r/7) ~ —1.802. Correspondingly, in the inter- 
val 1 < Q < 4 the locus Bq crosses the positive Q 
axis at the squares of these values, 4cos 2 (7r/(2^+ 1)) for 
£ = 0,1,2,3, i.e., 4, 1, 4cos 2 (tt/5) = (3 + VE)/2 ~ 2.618, 
and 4cos 2 (7r/7) ~ 3.247. For L x = 26, there are zeros 
very close to Q 4 = 2, Q 6 = 3, and Q$ = 2 + y/2 ~ 3.414. 
Our exact results also show that the locus B v contains 
a semi-infinite line segment on the real axis extending 
leftward from v = —2 to v = oo and, correspondingly, 
Bq contains a semi-infinite line segment on the interval 
4 < Q < oo. While the region boundary that includes 
Q = 1 and Q = 4 was convex in the neighborhood of 
Q = 1 for L y — 1 and L y — 2, it becomes concave in 
this neighborhood for L y = 3 (and L y = 4) . This bound- 
ary also extends further outward from the real axis than 
for the lower values of L y , reaching to roughly ±4.23i. 
The outermost finite complex-conjugate arcs reach far- 
ther away from the real axis than the analogous arcs for 
L y = 2, crossing the imaginary axis at Q ~ ±10i, and 
extend farther to the left, with arc endpoints at approx- 
imately Q ~ — 6.8 ± 4.8i. The dominant terms are as 
follows: (i) two different \3,i,max's in the region con- 
taining the real interval Q > 4 and the separate region 
extending to complex infinity, (ii) \3,2,max in the region 
containing the interval Q3 < Q < Q5, (iii) \3,3 t m ax in 
the region containing the interval Q5 < Q < Qj, i.e., 
approximately 2.618 < Q < 3.247, and (iv) ^4 = 1 in 
the region containing the interval Q7 < Q < Qi = 4. 
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FIG. 5: Locus Sq for the Potts model on a 3x oo graph with CDBC 



and Q 



This locus includes a semi-infinite line segment in the 



interval 4 < Q < oo. 
of length = 26. 



Partition function zeros are shown for a graph 



negative real v axis at v = — 2 and v = — 2 cos(7r/(2£+ 1)) 
for ^ = 1, 2, 3, 4 and thus Bq crosses the positive Q axis at 
the squares of these values. The outer complex-conjugate 
arcs on Bq continue the trend observed for L y = 2, 3 
of having parts that extend farther from the real axis, 
crossing the imaginary axis at Q ~ ±13i, but at the 
same time, having endpoints that move leftward and in- 
ward toward the real axis, so that these endpoints occur 
near to this axis. The dominant terms are as follows: 
(i) different A^i^a^'s in the region containing the in- 
terval Q > 4 and the region extending to complex in- 
finity, (ii) \/±.2,max in the region containing the interval 
Q3 < Q < <95) (hi) ^4,3,max in the region containing the 
interval Q5 < Q < Q7, (iv) X^A.max in the region con- 
taining the interval Q7 < Q < Qg, and (v) A^s, = 1 in 
the region containing the interval Qg < Q < Q\ = 4. In 
addition to the isolated zero at Q = Q2 = 0, the parti- 
tion function also exhibits other isolated zeros not on Bq 
very close to Q4 = 2, Qq = 3, Qs = 2 + \/2 ~ 3.414, and 
Q10 = (l/2)(5 + v / 5) — 3.618. The reason is again due to 
vanishing coefficients of dominant terms at these values 
of Q, as discussed above. 



V. GENERAL PROPERTIES AND 
CONJECTURE FOR B FOR ARBITRARY L v 




In our exact calculations on CDBC strips with widths 
1 < L y < 4, we have discerned several systematic fea- 
tures which we generalize here. For the Potts model on 
self-dual cyclic strips of the square lattice of arbitrary 
width L y , with Q and v satisfying eq. (|1.4|) . in the limit 
L x — > 00, we conjecture that 

• In the interval — 2 < v < — 1 the locus B v crosses 
the negative real v axis at v = — 2 and 



v = V 2 l+1 



-2 cos 



21 + 1 



for 1< I < L-n 



(5.1) 



For odd L y this locus B v also includes the semi- 
infinite line segment 

-oo<v<-2. (5.2) 



FIG. 6: Locus B plotted in the l/Q plane, for the Potts model on 
a 3 X 00 graph with CDBC and Q = v 2 . This locus includes a line 
segment in the interval < l/Q < 1/4. Partition function zeros 
are shown for a graph of length L x = 26. 



E. L v 



also 
,2 



calculated 



For the present work we have 
Z{Gcdbc[^*L x ], Q, v) exactly for Q = v 2 . Owing to the 
large number of A^dj's (126 in total), we again do not list 
them here. We find that the resultant locus B v crosses the 



• Hence, in the interval 1 < Q < 4 the locus Bq 
crosses the Q axis at 



Q = Q21+1 = 4 cos 



2^+1 



for < t < L„ 



(5.3) 



corresponding to the roots of unity given in eq. 
(f^TTTfl with r = 21+1. For odd L y the locus Bq 
also includes the semi-infinite line segment 



4 < Q < 00 



(5.4) 
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These properties do not preclude the possibility 
that, for sufficiently wide strips, the locus Bq may 
also cross the real Q axis at points in the intervals 
Q < or Q > 4. 

• For a strip of width L y , the locus B v divides the 
v plane into at least 2L y regions, and the locus 
Bq divides the Q plane into at least 2L y regions. 
(Note that the number of regions into which B v 
and Bq divide their respective planes need not be 
equal; for a gedanken example, if B v were to contain 
just two complex-conjugate arcs emanating from 
v = —2 and crossing the imaginary v axis at ±i&o 
with bo 7^ and then terminating in endpoints, this 
locus would not enclose a region in the v plane, but 
its image in the Q plane would enclose a region.) 

• In the respective regions of the v and Q planes 
away from the real axis and extending to complex 
infinity, the dominant eigenvalue is \L y ,i,max- In 
the region of the v plane containing the interval 
— oo < v < —2, and equivalently, in the region of 
the Q plane containing the interval 4 < Q < oo, 
the dominant eigenvalue is a different At, i, max . 
The fact that this is different is evident from the 
equimodular curve separating this region in the v 
plane from the region extending to complex infin- 
ity, and similarly in the Q plane. For odd L yi two 
eigenvalues from this set becomes degenerate on the 
real axis in the respective intervals — oo < v < — 2 
and 4 < Q < oo, giving rise to the line segments on 
B v and Bq. 

• Concerning the other regions that contain intervals 
of the real v and Q axes, the L y = 1 case has been 
solved exactly above: in the region — 2 < v < — I 
and the corresponding region 1 < Q < 4, = 1 
is dominant. For L y > 2, in the region of the v 
plane containing the interval V2k+3 < v < V2k+i, or 
equivalently in the region in the Q plane containing 
the interval Q2k+i < Q < Q2k+3 for 1 < k < L y — 1, 
Al ,k+i,max is dominant. In the respective regions 
of the v and Q planes containing the intervals 
-2 < v < v 2 L y +i and Q 2 l u +i < Q < Qi = 4, 
the dominant eigenvalue is \l v .l v +i = 1- 

• All of the Ai^^.j's have the magnitude |Aj, y ,d,j\ = 1 
at v — —2 or, in terms of Q, at Q = 4. Hence, all 
of the curves on B v meet at the point v — — 2, and 
correspondingly, all of the curves on Bq meet at 
the point Q = 4. 

• We next remark on zeros of Z(Gcdbc[L v x 
L x ],v 2 ,v) that, in the limit L x — > oo, do not lie 
on B, but instead are isolated (discrete). We have 
shown (so this is not just a conjecture) that there 
is an isolated zero at v — of multiplicity at least 
n + 1 and hence also at Q — with multiplicity 
that scales like n/2 for large n. In addition, from 
our exact results, we conjecture that for arbitrary 



L y , Z(G C DBc[L y x L x ],Q,v) with Q = v 2 has ze- 
ros close to the values 

Q = Q 2 e for 1 < I < L y + 1 . (5.5) 

The basis of this inference and the way to under- 
stand this result is that that these points are lo- 
cated within the region where the dominant eigen- 
value of the transfer matrix, \L y .d,max, has a coef- 
ficient in eq. I|3.2|) that vanishes at the points 
H5.5(l . as is evident from eqs. (|3.7|l and (|3.8|l . 

For each of our plots showing the exact loci Bq on 
infinite-length strips of the square lattice with cyclic self- 
dual boundary conditions and Q = v 2 , we have also 
shown zeros calculated for long finite strips for compari- 
son. As was evident, these zeros lie close to the asymp- 
totic loci B. This comparison is especially easy to make 
for the plots of B in the 1/Q plane, where it is compact; 
clearly, in the Q plane, the zeros for long finite sections 
cannot have arbitrarly large modulus and hence cannot 
track the portions of B that extend infinitely far from the 
origin. It is also of interest to compare our results with 
calculations of zeros of Z(Gcdbc[L x L], Q,v) for finite 
L x L sections of the square lattice with Q = v 2 in Ref. 
[l9|. Of course, for such a finite section, no locus Bq is 
defined; one is only able to make a rough comparison of 
patterns of zeros. One sees that a number of zeros in the 
Q plane occur at or near to certain Q r 's [19J and that the 
zeros in the v and Q planes exhibit patterns suggesting 
the importance of the points v = —2 and Q = 4. 

For general Q and v, the locus Bq is compact for the 
infinite-length, finite-width self-dual strips of the square 
lattice [U. We noted above how these strips share this 
property, along with their self-duality, with the (infinite) 
square lattice. In contrast, with these strips, when one 
restricts Q and v to satisfy the relation (|1.4|) , B v and Bq 
become noncompact in the respective v and Q planes, 
passing through 1/v = and 1/Q = 0. 

VI. SELF-DUAL STRIPS WITH FREE 
LONGITUDINAL BOUNDARY CONDITIONS 

We have also calculated Z(G,v 2 ,v) on self-dual strips 
of the square lattice with free longitudinal boundary con- 
ditions. A strip of this type of length L x and width L y 
is constructed by adding an external vertex which is con- 
nected to each of the vertices on one horizontal side, say 
the upper one, and to each of the vertices on one verti- 
cal side, say the right-hand one, of the strip [H[llli2- 
(Hence, the upper right-most vertex of the strip is con- 
nected to the external vertex with a double edge.) We 
denote this type of boundary condition as "free (self)- 
dual boundary condition", abbreviated as FDBC, and 
the graph as G FDB c[L y x L x \. 

For the strip Gfdbc\1 x m] we use our previous exact 
calculation of Z [2(J, |2l| and restrict Q and v to satisfy 
(II. 4|) . The partition function depends only on (the m'th 
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powers of) the \x,ij in cq. i|4.2[l . As is evident from eq. 
1)4. 2[1 . these are equal in magnitude for real v < —5/4, 
and consequently, B v is the semi-infinite line segment ex- 
tending from v — — 5/4 to u = — oo. The image of this in 
the Q plane, Bq, is the semi- infinite line segment from 
Q = (5/4) 2 = 1.5625 to Q = oo. Thus, as was the case 
with the infinite-length CDBC strips, the loci B v and Bq 
are noncompact in the respective v and Q planes. For 
FDBC strips with larger widths we find that B v and Bq 
consist of arcs and line segments and generically exhibit 
the above noncompactness. 



VII. OTHER STRIPS 

We have considered strips of the square lattice with 
free (F) and periodic (P) boundary conditions in the 
transverse and longitudinal directions, which we label 
in an obvious manner, e.g., (FBC y ,FBC x ) — free, 
(FBC y ,PBC x ) = cyclic, (PBC y ,FBC x ) = cylindrical, 
and (PBCy, PBC X ) = toroidal. We observe that for 
the cyclic square-lattice strips that we have studied, Bq 
crosses the real axis at a set of points including Q\ = 4 
and Q2£ for 1 < £ < L y . We find that for some strips Bq 
also crosses the real axis for Q < or Q > 4. 

In addition, we have calculated the loci B for infinite- 
length strips of the triangular lattice, with Q and v re- 
stricted to satisfy the relation Q = v 2 (v + 3), which is 
the phase transition condition for the Potts ferromagnetic 
on the triangular lattice. As in the case of the cyclic and 
toroidal square-lattice strips, we find that for the case 
of periodic boundary longitudinal conditions, Bq crosses 
the real axis at Q = and Q — 4 as well as at other points 
that depend on the type of strip. The intersections of Bq 
with the real axis sometimes include finite line segments. 
A detailed discussion of our results for these other lattice 
strips will be presented elsewhere. 



VIII. RELATION TO B FOR GENERAL Q, v 

It is interesting to observe the connections between the 
loci that we have determined for the Q = v 2 case to the 
continuous accumulation loci B in other variables, such 
as (i) complex v for fixed Q; (ii) complex Q for fixed v; 
and (iii) complex ( for, e.g., various values of Q. We 
concentrate on the CDBC strips and comment briefly on 
others. In Ref. 0| we displayed the loci in the £ plane 
for several values of Q in Ref. In relating these 

results to our present ones, one must take account of the 
noncommutativity i|3.11[l . Let us consider, for example, 
the CDBC strip of width L y = 1. Figs. 9, 11, 12, and 
13 of that paper gave Bq for Q = 2, 4, 5, and 100, calcu- 
lated by first setting Q equal to the given value, and then 
taking L x — ► oo; these plots show that this locus passes 
through f = — 1, and hence our special case Q — v 2 . 
This agrees with our present finding that for this special 
case, Bq passes through Q = Q\ = 4, and contains a 



semi- infinite line segment on the interval 4 < Q < oo. In 
the case Q = 2, the fact that the locus was calculated in 
Ref. [2l| by first setting Q equal to this value and then 
taking L x — > oo explains why the locus Bq in Ref. poj 
passes through £ = — 1, while our locus Bq, calculated 
by first taking L x — > oo, does not pass through Q = 2; 
the reason is that these limits do not commute. The ori- 
gin of this noncommutativity is clear: setting Q = 2 first 
means that = 0, and this annihilates what would be 
the dominant term in the vicinity of v = — V2, Q = 2, 
namely, (Ai :2 ) La! = 1- For Q — 3, the locus Bq, shown in 
Fig. 10 of Ref. [2lJ does not pass through ( = — 1, which 
again agrees with our finding that the point Q = 3 is not 
on Bq for the special case Q = v 2 . Similar comparisons 
can be made for the results shown in Ref. [2l| for CDBC 
strips with widths L y — 2 and L y = 3. 

The case of self-dual strips of the square lattice with 
free longitudinal boundary conditions was also studied 
in Ref. |2lJ. For example, as is shown in Figs. 2-6 of 
Ref. [IJ, as Q increases from 1 to (5/4) 2 , two complex- 
conjugate arcs on Bq enlongate and two of their endpoints 
come together and pinch the real axis at ( = —1; for 
Q > (5/4) 2 , the locus Bq continues to pass through ( = 
— 1. Again, this agrees with our present analysis, where 
we have found that for Q = v 2 , B v and Bq contain the 
respective semi-infinite line segments —00 < v < —5/4 
and (5/4) 2 < Q < 00. Similar connections can be made 
for other strips. 

We can also relate our present results to the loci Bq 
that we previously determined in the Q plane for the 
L x — > 00 limits of chromatic polynomials P(G, Q) = 
Z(G,Q,—1) of various lattice strip graphs with CDBC 
and FDBC [U- For CDBC strips we found that B 
crossed the real Q axis at Q = Q3 = 1, which lies on 
the submanifold Q — v 2 and agrees with our present 
result that for Q and v restricted to this submanifold, 
B v passes through v = — 1 (see Figs. 2-4 of Ref. 20] ). 
For strips of the square lattice with cyclic and toroidal 
(and Mobius and Klein-bottle) boundary conditions, we 
found that for v = — 1, Bq does not pass through Q = 1. 
This agrees with our finding that for these strips, with 
Q = v 2 , B v and Bq do not pass through the respective 
points v = — 1 and Q = 1. 



IX. CONCLUSIONS 

In conclusion, we have presented exact results for the 
continuous accumulation set B of the locus of zeros of 
the Potts model partition function for the infinite-length 
limits of self-dual cyclic square-lattice strips of various 
widths L y with Q and v satisfying the relation Q = v . 
For these quasi-one-dimensional strips, one motivation 
for interest in this submanifold of values of Q and v is 
its invariance under a duality transformation. From our 
exact calculations we have discerned several general fea- 
tures and have incorporated them in a conjecture appli- 
cable for arbitrary width L y . We have also presented 
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some results on B for self-dual strips of the square lat- NSC-94-2112-M-006-013 (S.-C.C). 
tice with free longitudinal boundary conditions. Further 
studies of these loci for other strips are worthwhile. The 
findings add to the set of exact results that one has for 
the Potts model on the n — > oo limits of lattice graphs. 
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